Breakdown of Scaling in the Nonequilibrium Critical Dynamics of the 

Two-Dimensional XY Model 
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The approach to equilibrium, from a nonequilibrium initial state, in a system at its critical point 
is usually described by a scaling theory with a single growing length scale, £(£) ~ t , where z is 
the dynamic exponent that governs the equilibrium dynamics. We show that, for the 2D XY model, 
the rate of approach to equilibrium depends on the initial condition. In particular, f (£) ~ t 1 ' 2 if no 
free vortices are present in the initial state, while £(£) ~ (t/lnt) ' if free vortices are present. 



While the theory of equilibrium critical phenomena has 
been a mature subject for more than 20 years, nonequilib- 
rium critical phenomena still pose some interesting chal- 
lenges. The simplest scenario consists of a system evolv- 
ing at its critical point from a nonequilibrium initial state 
in which the system was prepared at time t = 0. Since the 
characteristic relaxation time is infinite at criticality, an 
infinite system will never reach equilibrium. Instead, the 
system evolves towards equilibrium through a nonequilib- 
rium scaling state. Consider, for example, the equal-time 
pair correlation function, C(r,i) = (<f>(x., t) <f>(x + r,i)), 
where </> is the order-parameter field. In the nonequilib- 
rium scaling state it has the form 
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where d is dimension of space, r\ the usual critical expo- 
nent, and c is a constant. The scaling form (fil) holds in 
the limit r ^> a, £(£) ;§> a, with r/£(t) arbitrary, where a 
is a microscopic cut-off, e.g. a lattice spacing. The first 
factor in (Fy) is the equilibrium correlation function: Re- 
quiring that this be recovered for t = oo forces /(0) = 1. 
The physical interpretation of £(t) is the length scale 
up to which critical correlations have been established 
at time t: C(r,t) ~ cr -( d - 2 + r i) ( the equilibrium result 
holds, for a <C r <C £(£). Dynamical scaling suggests 



e(<) ~ t 1/z 



(2) 



for large t, where z is the usual dynamic exponent charac- 
terising temporal correlations in equilibrium. This result 
has been demonstrated in an expansion in e = 4 — d using 
standard field-theoretic renormalization group methods 
|ij. The importance of this result is that it shows that 
relaxation to equilibrium is governed by the same expo- 
nent as correlations in equilibrium. A second important 
result of reference fil is that the relation £(£) ~ t x l z holds 
independently of the nonequilibrium initial state, which 
can affect the scaling function, fix), in (nl) but not the 
exponent z (since this is a property of the equilibrium 
renormalization group fixed point). 

Two special cases illustrate the dependence of f(x) on 
the initial conditions. For a disordered initial condition, 
the system will remain disordered on scales r 3> £(£), 



so f(x) will fall off rapidly for x 3> 1. For an ini- 
tial condition with long-range order (i.e. non-zero initial 
magnetization) , dynamical scaling predicts that the mag- 
netization M(t) will decay asymptotically as t~ l3 ^ l/z ~ 
£,{t)~P' v - In this case C{r,t) approaches C(oo,i) = 
M 2 (t) - t~ 20/vz = t-^ d ' 2+ri ^ z using standard scaling 
laws. So in this case f(x) ~ x d ~ 2+r > for x — > oo. 

The purpose of this Letter is to challenge this sim- 
ple picture for the XY model in d = 2 at (and below) 
the Kosterlitz-Thouless (KT) transition 0]. Specifically 
we argue that the growing length scale £(£) satisfies (0), 
with z = 2, if the initial state contains no free vortices, 
whereas £(£) ~ (i/lni) 1 / 2 if free vortices are present. We 
first present numerical simulation results supporting this 
scenario, and then provide the theoretical interpretation. 
The two types of initial state we shall consider are (a) 
completely ordered (no free vortices) , and (b) completely 
random (free vortices present). 

The XY model consists of planar spins {Si} at the 
sites of a square lattice of linear size L, with Hamilto- 
nian H = — J2<% ? > >% ' Sj> where the sum is over lattice 
links and we have taken the exchange interaction to have 
strength unity. We adopt conventional 'heat-bath' dy- 
namics in which a spin is moved to a trial configuration 
chosen at random on the unit circle, and the move ac- 
cepted with probability [1 + exp(A£;/T)] _1 , where AE 
is the energy change associated with the move, and T is 
the temperature. The lattice is divided into two sublat- 
tices, and the sublattices updated alternately. One unit 
of time corresponds to an attempted move of every spin. 

A convenient quantity to study is the 'time-dependent 
Binder cumulant' j|,|J, <7i,(t), defined by 
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where M(t) — ^\ Si(t) is the total magnetization at time 
t, and (. . .) indicates an average over independent Monte 
Carlo runs (I0 4 runs were used in practice). Because 
the powers of M in numerator and denominator are the 
same, gh{t) depends (at a critical point, and if dynam- 
ical scaling holds) only on the ratio £(t)/L, where L is 
the (linear) size of the lattice: 
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provided, as always, that both £(£) and L are suffi- 
ciently large. This result provides the basis for a de- 
termination of £(£) using finite-size scaling. If conven- 
tional dynamical scaling holds, the scaling function G(x) 
may depend on the initial state, but the scaling variable, 
£,{€)/ L = t x l z I L, will not. For an ordered initial state (all 
spins parallel) G(0) = 1, while for a random initial condi- 
tion (each spin chosen independently from a unit circle) 
G(0) = follows from the Gaussian distribution (central 
limit theorem) of M(0). For t — > oo, G(x) approaches, 
in both cases, the universal value G(oo) characteristic of 
the critical point. For the KT phase there is actually a 
line of such fixed points, T < Tkt, (and a corresponding 
set of values Gt (oo)) but we will focus primarily on the 
KT point, Tkt, using the accepted value Tkt = 0.90 |y. 
Data for the ordered initial state are presented in Fig- 
ure 1, for system sizes L = 12, 16, 24, 32, and 48. The ab- 
scissa, t/L 2 , corresponds to a scaling variable t/L z with 
z = 2. This choice of z is dictated by the spin- wave the- 
ory (i.e. no free vortices) that describes the large-scale 
properties of the KT phase everywhere along the fixed 
line T < Tkt- The best collapse using all the data fa- 
vors a slightly lower value, but the value z — 2 clearly 
gives a good scaling collapse for larger L, i.e. L > 24 
(note the expanded scale compared to Figures 2 and 3). 
Collapsing the data for pairs of L gives effective expo- 
nents z(Li,L 2 ) given by z(12, 16) = 1.75(5), z(16,24) = 
1.83(3), z(24,32) = 1.96(2), z(32,48) = 2.00(2), consis- 
tent with a convergence to z — 2 for L — ► oo. Recent sim- 
ulations by Luo et al. give similar results: z — 1.96(4) 
for T — 0.90 and an ordered initial state. 
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FIG. 1. Scaling plot, with z — 2, for the time-dependent 
Binder parameter, starting from an ordered initial condition, 
for system sizes L = 12, 16, 24, 32, 48 



The date for a random initial condition are presented in 
Figures 2 and 3. In Figure 2, we attempt to collapse the 
data with a scaling variable t/L z . The scaling collapse 
is very good, but a much higher value of the dynamical 
exponent, z ~ 2.35, is required than for an ordered initial 
state. For a random initial condition Luo et al. found, 
by direct measurement of the time-dependence for a large 
lattice (L — 512), the slightly smaller result z — 2.29(1) 




FIG. 2. Scaling plot, with z — 2.35, for the time-dependent 
Binder parameter, starting from a disordered initial condition, 
for system sizes L = 12, 16, 24, 32, 48 



At first sight, these results seem remarkable: Different 
values of z are required to fit the approach to equilib- 
rium from ordered (or 'low-temperature') and disordered 
(or 'high-temperature') initial states, whereas dynamical 
scaling predicts a unique value of z, namely that which 
describes equilibrium correlations (in this case z = 2). 
What is going on here? It is worth noting that for Ising 
systems, the two different initial conditions give compat- 
ible results J7|. The data for the XY model in d = 2 
seem to point clearly to a breakdown of dynamical scal- 
ing. This is indeed our conclusion, but the breakdown 
is weaker than the naive fit shown in Figure 2 suggests. 
We will argue that, for a disordered initial condition, the 
characteristic length scale £(£) grows as (£/ lni) 1 ' 2 rather 
than t x / z . Before presenting the arguments, we test this 
prediction in Figure 3, where t/L 2 ln(t/io) is used as ab- 
scissa. The fit is excellent. The value to = 0.5 was used 
for the short-time cut-off, but the fit is not too sensitive 
to this value. 

The quality of the scaling collapses in Figures 2 and 
3 are comparable, but the fit used in Figure 3 has a 
theoretical underpinning. First, however, we note that 
the scaling form (pi), with z = 2, follows from the spin- 



wave theory for an ordered initial state: No free vor- 
tices are present at t — 0, and none gets generated by 
thermal noise for any T < Tkt- The calculation of 
C(r, t) is straightforward [p|, and gives precisely the scal- 
ing form (|) with £(t) = i 1 / 2 and f{x) = exp[-7j J(x)/2], 

where J(x) — f^ (dy/y) [1 — exp(— y)] and 77 = 1/4 for 
T = Tkt- 
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FIG. 3. Same as Figure 2, but for scaling variable 
i/[L 2 ln(t/t )], with to = 0.5. 



For a disordered initial condition, very different con- 
siderations are involved. The initial state contains many 
free vortices and antivortices. The approach to the equi- 
librium critical state proceeds through the annihilation 
of vortex-antivortex pairs, which is a slower process than 
the equilibration of spin waves. For pedagogical pur- 
poses, we consider first the case where the system evolves 
at T = 0, instead of Tkt- The evolution of the sys- 
tem via vortex-antivortex annihilation is an example of 
phase-ordering dynamics |9j. It is convenient to adopt 
a continuum approach based on the non-linear sigma 
model Hamiltonian H = (1/2) J d 2 r (V0) 2 , with local 
constraint <fi 2 — 1 . A field configuration describing a sin- 
gle free vortex, (j) = r/\r\, has an energy E v — n\ii(L/a), 
where L and a are the system size and microscopic cut- 
off as before. A vortex-antivortex pair, separated by dis- 
tance R, screen each other's far fields at scales larger 
than R, leading to a pair energy E p ~ 27rln(_R/a), and 
an attractive force F = —dE p /dR = —2tt/R between the 
vortex and the antivortex. 

To discuss pair annihilation, some dynamics has to 
be imposed. The Monte-Carlo dynamics used here is 
in the 'non-conserved' universality class (i.e. the magne- 
tization is not conserved) described (at T = 0) by the 
continuum model d<p/ dt — —SH/S(f>- This equation can 
be used |10| to compute an effective friction constant 
7?(i?) associated with the motion of the vortex and an- 



tivortex under the force F . An isolated vortex mov- 
ing at speed v in the ^-direction has field configuration 
4>{x,y,t) = 4> v {x — vt,y). Energy is dissipated at a rate 
dE/dt = / d 2 r (5E/6$) ■ {d$/dt) = - J d 2 r (d$/dt) 2 = 
—v 2 J d 2 r (d(j) v j dx) 2 = —r/ v v 2 . Inserting the equilibrium 
vortex configuration, which is isotropic, gives the limit- 
ing zero- velocity friction constant as 770 = E v , i.e. 770, like 
the vortex energy E v , diverges logarithmically with the 
system size, L. For a vortex-antivortex pair, this trans- 
lates into a logarithmic dependence on the separation 
@,77(i?) = 7rln(i?/a). 

In the many-vortex situation envisaged for the 
nonequilibrium critical dynamics, the usual scaling ar- 
guments [pr[y] can be invoked, in which the pair sep- 
aration, R, is replaced by the typical spacing, £(t), be- 
tween vortices and antivortices. The typical force on a 
vortex (or antivortex) is then F ~ l/£, while the typ- 
ical friction constant is r\ ~ ln(£/a). so the typical 
speed of a vortex is d^/dt ~ Fjr\ ~ l/[£ln(£/a)], giv- 
ing £(£) ~ [t/ln^t/to)] 1 ' 2 , with t ~ a 2 . An alternative 
approach leading to the same result is given in [[Ll| . 

For all T in the range < T < Tkt, the large- 
scale properties in equilibrium are controlled by a fixed 
point with zero vortex fugacity, i.e. by the spin- wave the- 
ory, where the role of bound vortex-antivortex pairs is 
to renormalize the spin- wave stiffness. In the nonequi- 
librium case where free vortices and antivortices are 
present, due to a disordered initial condition, the dy- 
namics on scales less than £(£) should therefore be de- 
scribed by renormalized spin-wave theory, and the result 
£(i) ~ [t/ln(t/to)] ly ' 2 should apply to all temperatures 
T < Tkt, including Tkt itself. This is our interpreta- 
tion of the data in Figure 3. It accounts for the good 
data collapse using the appropriate scaling variable. 

Clearly our result, £(£) ~ (t/lni) 1 / 2 for a disordered 
initial condition, is asymptotically equivalent to an expo- 
nent z = 2 (though the logarithmic correction still repre- 
sents a scaling violation). So as L and t are increased we 
would expect the effective exponent, obtained by forcing 
a fit with a scaling variable t/L z , to decrease towards 
2. Collapsing data for pairs of L gives effective expo- 
nents z{L ll L 2 ) given by z(12, 16) = 2.47(3), z(16,24) = 
2.37(3), z(24, 32) = 2.29(3), and 2(32, 48) = 2.34(3). The 
quoted errors are subjective. They are estimated from 
the quality of the data collapse, but make no allowance 
for statistical errors in the data. They therefore represent 
lower bounds on the true errors Q . With this caveat the 
overall decreasing trend of the effective z with increasing 
L is clear, and accords with our expectations. 

Although the data presented here are restricted to 
the Kosterlitz-Thouless transition temperature, Tkt, the 
theoretical interpretation we have outlined holds for all 
T < Tkt- In reference H, data were obtained for a range 
of temperatures at and below T KT : T = 0.90, 0.86, 0.80, 
and 0.70. For a uniform initial state, the corresponding 



values of z are 1.96(4), 1.98(4), 1.94(2), and 1.98(4), con- 
sistent with the result z = 2, for all T < Tkt, expected 
from spin-wave theory. The equivalent effective expo- 
nents obtained with a disordered initial state are consis- 
tently larger: z = 2.29(1), 2.31(2), 2.33(1), and 2.38(2). 
We have argued that the correct interpretation of these 
anomalously large z values is a logarithmically modified 
growth, £(£) ~ [t/ln^/io)] 1 ' 2 ) of the characteristic length 
scale. The slow increase of the effective exponent with 
decreasing T can be accounted for by a weak temperature 
dependence of the time-scale to inside the logarithm. 

For the ordered initial condition, the scaling function 
G(x) in equation (4) can, in principle, be calculated ex- 
actly using the spin-wave theory. This is technically more 
difficult, however, than the calculation || of the pair cor- 
relation function, because the evaluation of ((Af 2 ) 2 ) in- 
volves 4-point correlation functions. We hope to present 
a detailed theory for G{x) in future work. 

To summarize, we have argued that the rate of ap- 
proach to equilibrium at (and below) the Kosterlitz- 
Thoulcss transition temperature depends on whether or 
not the initial state contains unbound vortices. Thus for 
a disordered initial state, where free vortices are present, 
the relaxation to equilibrium is slower, by logarithmic 
factors, than for an ordered initial state where no free 
vortices are present. It is possible that this result is pe- 
culiar to systems with defect-driven phase transitions. It 
goes against the expectation [Q that the scale length £(£) 
controlling the relaxation to equilibrium, e.g. in (1) and 
(4) , should be independent of the initial conditions (al- 
though the corresponding scaling functions may not be). 
This expectation is based on a perturbative renormal- 
ization group treatment in 4 — e dimensions. Such an 
approach is not, presumably, sensitive to the effects of 
topological defects (vortices, in this case), which are the 
source of the scaling violations reported here. 

A. Bray thanks M. A. Moore for a helpful discussion. 
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